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CN ■ Abstract 

We investigate femtoscopic tttt correlations using the UrQMD approach com- 
bined with a correlation afterburner. The dependence of vrvr correlations on 
Q ^ the charged particle multiplicity and formation time in p+p collisions at 

^ ■ \/snn= 7 TeV is explored and compared to present ALICE data. The data 

allows to constrain the formation time in the string fragmentation to r/ < 0.8 
fm/c. 

^ ! 1. Introduction 

'NT . 

• ! With the start of the LHC physics program two years ago a tremen- 

Q I dous amount of new data became available. Apart from the heavy ion data, 

CN I the proton-proton (pp) program allowed to explore collective features of the 

strong interaction in high multiplicity pp events. 

It seems like in massive nucleus-nucleus collisions, a strongly interacting 
medium is created even in pp collisions, that exhibits similar bulkproperties 
j^ I such as space momentum correlations and collective behaviour [ll, |2|, y, y]. 

The details of these correlations can be probed using Hanbury Brown-Twiss 
(HBT) |5| interferometry. While it is often argued, that the particle emitting 
system in p-t-p collisions is too small to create a medium that exhibits bulk 
properties, this is different at a center of mass energy of a/s= 7 TeV [6[ . Here, 
the particle multiplicity is about the same as in nucleus-nucleus collisions, 
studied at the Relativistic Heavy Ion Collider (RHIC) in Brookhaven. For 
previous studies of femtoscopic correlations in p+p collisions at RHIC and 
Tevatron see [Tl, |8|. This data suggests that space momentum correlations 
are developed even in pp collisions as soon, as high particle multiplicities are 
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achieved. Thus, it is worthwhile studying the dependence of HBT observables 
on the event multiphcity. As the system created in p+p coUisions at LHC 
is still small, an essential quantity that influences the particle freezeout radii 
is the formation time in flux tube fragmentation. Without going into the 
details of the specific implementation, it is clear that the formation time sets 
the scale for a minimum value of the source lifetime - of course followed by 
resonance decay and rescattering. In this paper we suggest, that the recent 
LHC data on pp collisions allows to determine the formation time in the flux 
tube break-up. Results for Pb+Pb reactions and scaling studies at the LHC 
within the same model can be found in J9l. llO . 



2. Two-particle correlations 

To this aim we apply two-particle correlations to extract the space-time 
structure of ultra- relativistic particle collisions. For bosons the basic equation 
for the correlation function C(q) is 

Cq,K =1 + ^^^ ,2' , 1 

where q = pi — p2, K = i(pi -|- P2), with pi, p2 being the four-momenta 
of the particles and S(x,K) is the phase-space density of the particle emit- 
ting source. Using the Pratt-Bertsch parametrization we fit the resulting 
correlation function with the functional form 



C(q,K) = l + X{K)exp 



- J2 m^U^) 

. i,j=o,s,l 



(2) 



using the OSL-sjsteTa. Here, / is the longitudinal direction along the beam 
axis, the out direction along the transverse component of K and s the 
side direction perpendicular to the afore mentioned directions. Rfj(K.) are 
the HBT radii. All the radii in this paper are computed in the longitudinal 
comoving system (Icms). 

The width of the pion pair separation in space is inversely proportional 
to the momentum difference q probing it. In the presence of flow, high p^ 
particles that are far apart in position space are very likely to be also far 
apart in momentum space. This results in smaller HBT radii than the actual 



source size, because only the region of homogeneity is probed and not the 
entire source. For particles with low p^ the effect is not as strong, resulting 
in a larger probed region, thus giving larger HBT radii. 

3. Particle formation time 

The formation time denotes the time it takes for a hadron to be produced 
from a fragmenting string. A very common model to describe such a flux 



tube fragmentation is the Lund string model [11|. In the Lund model the 
formation time consists of the time it takes to produce an quark antiquark 
pair and the time it takes for that pair to form a hadron. For the Lund 
model both of these times are proportional to the transverse mass of the cre- 
ated hadron and inversely proportional to the string tension. For simplicity 
UrQMD uses a constant formation time of Tf = 0.8 fm/c for hard collisions. 
Only after the formation time particles, e.g. p-mesons, can decay or perform 
subsequent scatterings with their full cross section. During their formation 
time leading hadrons are allowed to interact with a reduced cross section. 
Other particles are not allowed to interact at all during their formation time. 
While a Heaviside function like behaviour of the cross section is implemented 



in UrQMD, there are other model studies |12l . |13] that investigate the influ- 
ence of several scenarios of continuous changes from zero to full cross section 
during the formation time. The freezeout space-time positions of hadrons 
are defined as either their point of formation or the point of last interaction, 
whatever occurs later in time. Since HBT probes the freezeout distribution 
the extracted radii are sensitive on the value of Tf in small systems at high 
collision energies where the majority of particles is produced from flux tube 
fragmentation. Although there are many theoretical studies on the formation 



time [IJ, |15|, |16|, ll7|, ll8| few of them allow to put constraints on Tf or the 



behaviour of the cross section growth from experimental data. 

4. Analysis process 

The UrQMD (v3.3pl) is used to generate the freezeout distributions of 
particles in proton-proton collisions at ^/sivm= 7 TeV. For further informa- 
tion about UrQMD the reader is referred to 19|, |20|, l2l| . After the simulation 



run, a correlation afterburner is applied to the freezeout distribution to calcu- 
late the corresponding 3D correlation function using Eq. [TJ In this analysis, 
as in ALICE data, particles within a pseudorapidity interval of \ri\ < 1.2 



charged 


UrQMD {dN,h/dr^)\^\^,,2 


ALICE {dN,h/dr^)\,\^,,, 


1-11 


2.52 


3.2 


12-16 


5.74 


7.4 


17-22 


8.01 


10.4 


23-29 


10.72 


13.6 


30-36 


13.65 


17.1 


37-44 


16.74 


20.2 


45-57 


20.94 


24.2 


58-149 


27.57 


31.1 



Table 1: Table of the investigated multiplicity intervals. The first column shows the 
interval boundaries, the second column the mean charged particle multiplicity per unit of 
pseudorapidity [dNch/dr]) from events with at least one charged particle in \ri\ < 1.2 from 
UrQMD. The third column shows the same quantity from ALICE data '22\. 



are taken into account. The analysis is done differentially for K± bins of 
100 MeV in the region of K± = 0.1 — 1.0 GeV and also for the different 
event multiplicities listed in Table [H Generally the average dNch/drj from 
UrQMD is 15% smaller than the one measured by the ALICE collaboration, 
in the same charged multiplicity classes, because we did not employ specific 
PYTHIA tunes for the present analysis. All the correlation functions are 
computed for pairs in the longitudinal comoving system (Icms). The HBT 
radii are extracted by fitting Eq. [2] to the 3-dimensional correlation functions 
over a range of \qi\ < 800 MeV. 

5. Results 



In this section the results on HBT radii from the UrQMD model are com- 



pared to ALICE data [22|. In Fig. [T]the projections in out, side and long 
direction of the 3D correlation function together with a projection of the 
fit for i^^=0.3-0.4 GeV and N'^^^ ,= 23-29 are shown as an example. It 
can be seen, that the calculated correlation function (shown as dots) is well 
described by a Gaussian fit (lines). However, oscillations of the correlation 
function are present at larger q. This indicates that there is a non-gaussian 
component in the underlying separation distribution of the pion freezeout 
points. Here we will not investigate this interesting question further. 
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Figure 1: The dots represent projections of the 3 dimensional correlation function for 
KjL^0.3-0A GeV and iv]y^Ltged= 23-29. The lines represent a x^ St of Eq. H to the 
correlation function. Both the result of the fit and the correlation function are integrated 
over a range of qi = ±0.17 GeV in the other directions for the purpose of projection. 



The K±^ dependence of the HBT radii extracted from the UrQMD freeze- 
out distribution is presented in Fig. |2]for the dNch/drj classes defined in Tab. 
[T]in the pseudorapidity interval |?7| < 1.2 in comparison to the ALICE data. 
The UrQMD calculations are presented for different values of the formation 
time Tf {Tf= 0.3 fm/c, dashed lines; rj= 0.8 fm/c, full lines; t/=2 fm/c, 
dotted lines). For Tf =0.3 fm/c one obtains a good description for Rout, while 
Rside is slightly under predicted and the values for Riong are in line with data 
from ALICE. The choice Tf =0.8 fm/c leads to a slight overestimation of Rout, 
however it leads to a reasonable description of Rside data. Also the K± be- 
haviour in Riong is much closer to the behaviour of the data. In contrast, a 
formation time of Tf=2 fm/c, leads to a drastic over estimation of the data 
for all radii. Although there are discrepancies between model and data for 
all values of Tf, the sensitivity on r/ is much larger than those discrepancies. 
Therefore, the present ALICE data allows to constrain the formation time 
to values of Tf ~ 0.3-0.8 fm/c. 
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Figure 2: The lines represent HBT radii in pp collisions at \/s^ 7 TcV from UrQMD 
for different multiplicities and formation times. The various line styles refer to results for 
Tf = 0.3 fm/c (dashed), Tf = 0.8 fm/c (default - full lines) and Tf = 2 fm/c (dotted). The 
colors represent the multiplicity classes as defined in table [1] The Points are data from 
the ALICE experiment |22| . 



Let US finally discuss the overall shape of the radii as a function of multi- 
plicity and K±. The Rside radii (see Fig. [21 middle) from UrQMD and in the 
data are fiat as a function of K± for low multiplicity events. With increasing 
multiplicity the radii develop a decrease towards higher K±. This is exactly 
the behavior one would expect for the development of space-momentum cor- 
relations with rising event multiplicity |6|. For Rout (Fig. HI left) however, 
there is a K± dependence present in all multiplicity classes. Thus, only the 
development of radial flow with rising particle multiplicity seems not suffi- 
cient to explain the K± dependence. Since Rout and Riong contain lifetime 
contributions of the source and Rside does not, there seems to be an addi- 
tional non-trivial K± and multiplicity dependence in the emission duration 
needed to explain the difference in the behaviour of Rout, Rside and Riong- 



This additional correlation might be due to an additional momentum depen- 
dence in Tf apart from the trivial Lorentz boost. This would lead to a direct 
effect on the emission duration, because it changes the particles production 
spacetime points. It would also lead to an indirect change of the emission 
region, since the particle rescattering is influenced by Tf. In this case of 
pp collisions the effect of the rescattering should be negligible. Preformed 



hadron interactions become important important in AA collisions [15|, |23| . 



6. Summary 

We have performed an analysis of proton-proton collisions at the top LHC 
energy of ^/s = 7 TeV. We have shown that the data provides rather direct 
access to the particle formation times in the flux tube fragmentation. The 
sensitivity to Tf is large enough compared to the model uncertainties to find 
that a value Tf = 0.3 — 0.8 fm/c is strongly favored compared to larger values 
for Tf. Values of rj > 2 fm/c can be ruled out from the present analysis. 
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